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We study theoretically the hydrodynamics of a fluid drop containing oriented filaments 
endowed with active contractile or extensile stresses and placed on a solid surface. The 
active stresses alter qualitatively the wetting properties of the drop, leading to new 
spreading laws and novel static drop shapes. Candidate systems for testing our predictions 
include cytoskeletal extracts with motors and ATP, suspensions of bacteria or pulsatile 
cells, or fluids laden with artificial self-propelled colloids. 



1. Introduction and results 

1.1. Background 



fluid mechanics ( Pedley & Kessler 


1992; Ramachandran et all 


20061 ICisneros et al. 2007; 


Mehandia & Nott 12008; Ishikawa 


20091; Lauea & Powers! 20091 iHernandez-Ortiz et al. 


20091 Koch & Subramanianl 2011 


). One line of recent progress in this field has come 



through the observation that ordered phases of active particles could be viewed as liquid 
crystals with a key novel feature: the constituent particles are endow e d with self-generated 
"active stresses" (iTqner et all l2005t iLiverpool fc Marchettil 120061; Ijiilicher et all 12007 



IJoannv fc Prostll2009llBaskaran fc Marchettill2009l : lRamaswamvll2010l : lMenonll2oio| ). The 
resulting extension of liquid crystal hydrodynamics to include self-propelling stresses. 



initia lly proposed as a coarse-grained description of swimmers (jSimha fc R amaswamv 



2002) , also emerged naturally in a theory of the cytoskeleton as an active gel (jKruse et al 



20041) . Within liquid-crystal hydrodynamics one could imagine a wide variety of orien- 
tationally or translationally ordered states, but most recent work has focused on the 
case of local uniaxial orientational order. In such situations it is useful to distinguish 
active fluids as contractile or extensile, according to whether the active stress gener- 
ates flow inward or outward along the local ordering axis. This classification corresponds 
to the puller-pu sher classification at the scale of a single swimming particle (see, e.g., 
Saintillanfe oiO). However, for systems such as motor-filament extracts, where a pulling 



or pushing motif is not obvious at the scale of a single particle, the sign of the ac- 
tive stress at a coarse-grained scale is best established by an independent measurement 
(jTanaka-Takiguchi et all\2004 ). A further distinction is necessary, between polar states 
which break fore-aft symmetry along the axis of orientation, and apolar states which 
are fore-aft symmetric. Swimmers in a macroscopically coherent flock are in a state of 
polar order. However, to leading order in inverse powers of distance from the particle, 
i.e., to force-dipole order, the flow-field generated by a general swimming particle is the 
same as that produced by one that merely stirs its surroundings without propelling itself. 
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The distinction between "movers" and "shakers" (jHatwalne et all\2004 ) lies in their near 
field and is thus subdominant in a gradient expans i on. When studying; rheology (see, e.g., 
Hatwalne et a/.ll2004 iLiverpool fe MarchettfeOOd; ICates et aI1l2008llSokolov fc Aranson 



2O10HRafai et aLll201o[Tsaintillanll2oioli " and large-scale flow, it is therefore a reasonable 
starting point to ignore the consequences of polar order in active fluids. Accordingly 
our study in this paper is restricted to states with apol ar uniaxial orienta tional order. 
The effect of polarity on the instability of active fluids ( Giomi et al. 20081 ). the conse- 
quences of bottom-h eaviness and more generally the interplay of swimming with gravity 
(jHill fe Pedlevll2005l ) all lie outside the scope of this work. 

Instability and self-generated flow are conspicuous features of active fluids. Indeed, the 
idea of hydrodynamic inst ability in living matter, and even the term active stress, goes 
back at least to the work of iFinlavson fc Scriven ( 19691 ) . Their primary interest, however, 
lay in cases in which the stress was built from tensorial combinations of gradients of 
local scalars such as concentration and temperature. The present work, and much recent 
work on flows in active matter, relies on the ordering of orientational degrees of freedom 
to transmit the effects of local activity to large scales. A bulk sample of active fluid 
or suspension in a state of uniaxial orientational order is under macroscopic uniaxial 
compression or tension, and must therefore buckle along or transverse to the axis. This 
is why the non-flow ing state of a bulk active fluid, whether extensile or contractile, is 
in general unstable (|Simha fc RamaswamvM2002t IVoituriez et aLl I2005T ) towards director 
distortion and flow without the assistance of externally imposed pressure gradients. 

Theoretical studies of active fluids in a thin-film geometry with one or more free sur- 
faces predict particularly stri king effects for a refe rence state in which the orientation 
field lies parallel to the film. Voituriez et al. ( 20061 ) consider a freely-suspended apolar 



film, ignore the displacements o f the free surfaces, and predict the sp ontaneous genera- 
tion of topological defect arrays. ISankararaman fc Ramaswamvl ( 20091 ) consider an active 



film spread on a solid surface, and show that the coupling of polar orientational order to 
the tilt of the free surface of the film leads to growing and propagating undulations. Both 
these studies consider films unbounded in the lateral plane. Finite drops of s uspensions 



of self-prope lled organisms have been studied in theory and experiment by iBees et 



(|2000ll2002|) . but the key ingredients in those studies were cell division, the secretion of 
a wetting agent, and concentration dependence of the viscosity. 

In this paper we examine theoretically the effect of self-propelling stresses on the shape 
and the spreading behaviour of a finite drop of active fluid on a solid substrate, under 
conditions of partial wetting with small equilibrium contact angle, and local polarization 
anchored parallel to the surface of the drop. This boundary condition is in contrast to 
the perpendicular alignment f requently encountered in the liquid-crystal literature; see, 
e.g.. lGupta fc Abbott! (|l999l ): lYamaguchi fc Sato! (|l996l k IPeGennes fc Prostl (| 19931 ) We 
consider separately the cases of contractile and extensile active stresses, and study two 
situations: steady standing drops, with vanishing velocity field, and spreading drops. 

1.2. Summary of results 

Here are our main results in brief, (i) The shape of the standing drop is determined pri- 
marily by the interplay of active stresses o~q and surface tension 7, with the liquid-crystal 
elasticity of the ordered filaments playing a role near the contact line. The stationary 
height profile h(x) as a function of the horizontal coordinate x can be obtained as the 
Euler-Lagrange equation for position h, as a function of "time" x, of a particle in an 
effective potential V(h) = —ao(hlnh/ho — h). (ii) For a 3-dimensional drop, the anchor- 
ing conditions impose the existence of defects in the polarization field. We propose three 
possible defect structures; both the final state of spreading and the spreading kinetics 
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Figure 1. A drop characterized by height and orientation fields h and p, the latter constrained 
to lie parallel to the bounding surfaces. For the present apolar case, p and — p are equivalent, 
so we depict the orientation by headless arrows. For clarity we have shown the filaments and 
the orientation field only near the edge of the drop. 




Figure 2. Sketches of the three defect configurations of the nematic director field that we 
consider: (a) the aster, (b) the vortex and (c) the pair of half-defects. 



depend on the precise nature of these defects, (iii) Extensile active stresses lead to flat 
drops, with thickness varying as 7# 2 /o"o where 9 is the equilibrium contact angle, (iv) 
The interplay between active and viscous stresses leads to anomalously rapid spreading 
of the drop as a function of time t. In the simple case where the orientation field has 
an aster-like defect line at the center of the drop, the drop retains rotational symme- 
try with diameter R ~ (aofl 2 t/ /i) 1 / 6 where fi and /i are its volume and shear viscosity 
respectively. For a 1 + 1-dimensional geometry, we predict a i 1 / 4 growth law. (v) For a 
vortex defect, in which the orientation axis points in the azimuthal direction, spreading is 
arrested at long enough times, and the final height profile of the drop is a non-monotone 
function of radius as shown in Eq. (|2.17[) and Fig. [5l For contractile active stresses the 
drop is fatter than in the absence of activity (Fig. HJ). 
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2. Detailed theory 

We now describe our calculations in more detail. Consider as in Fig. Q] a drop with a 
shape given by its thickness h(r±, t) in the vertical (z) direction as a function of in-plane 
position rj_ = (x, y) at time t and containing filaments with orient ation field p. The active 



character of the fluid en dows the drop with an active contribution (jSimha &: Ramaswamy 
2002; iKruse et al. I l2004h to the stress tensor of the medium. 



a a = <j pp (2.1) 

Negative and positive values of o~a correspond respectively to contractile and extensile 
active stresses, and the magnitude of Co is proportional to the local number density 
of filaments c, which for simplicity we assume to be constant and uniform, although 
interesting physical effects are likely to arise if the dynamics of c is taken into account. 
The anchoring condition is such that there is no normal component of the polarization at 
any bounding surface with which the suspension is in contac t. This boundary condition 
is motivated in part by observations (jKemkemer et a.71 |2000T) on thin films of amoeboid 
cells, in which the cells lie in the plane of the glass slide on which they are spread, and 
form nematic liquid crystalline structures. Although we use a vectorial orientation field, 
in this paper we only discuss apolar fluids where p and — p are equivalent. This arises 
if the filaments lack a head-tail distinction, or are ordered with axes parallel on average 
but without respecting the head-tail polarity. Our aim is to obtain equations of motion 
for h and a z-averaged orientation p(rj_,i), whose static solutions yield the steady-state 
shape of t he drop and whose dynamic scaling pro perties lead to the spreading laws for 
the drop (iLeeer fc Joannvlll992l: DeGenneslll985l ). Our treatment is along the lines of 



Bcn Amar fc Cummingsl (|2001l ). taking into account, as in lSankararaman fc Ramaswamy 



(|2009l ) , the presence of active stresses but ignoring the effects of polarity, and applied to 
a drop rather than an unbounded film. 

We solve for the hydrodynamic velocity field u = (uj_, u z ) inside the drop for a given 
thickness profile, and o btain the equations of motion for h from the kinematic boundary 
condition ( Stonei 20051 ) h — u z — u± ■ V±h connecting the thickness h to the velocity field 



u evaluated at the free surface. The condition of incompressibility V • u = turns the 
height dynamics into a local conservation law 

d t h + V ± • \ u ± dz = (2.2) 
Jo 

in the _L plane. For simplicity, we only consider here the case where the equilibrium 
contact angle of the drop 9 is small. As appropriate for slow, small-scale thin-film flows 
we obtain the thickness-averaged horizontal velocity in (|2.2[) in terms of p , c and h by 
solving the steady Stokes equation which, in the lubrication approximation ( Stone 2005) 
u z = 0, |Vj_u| <C |9 z u|, reads 

M<9fu_L - VP + V • (o- a + cr e ) = (2.3) 

where [i and P are the viscosity and pressure field of the fluid, a e is the contribution to 
the stress from the nem atic elasticity of p (see, e.g., Leslie! 19791 Chandrasekhar 1992 : 



DeGennes &: Prost 1993f l and a a is the active stress introduced above. 



2.1. 1+1- dimensional treatment 

It is instructive to solve the problem in a simple (l + l)-dimensional picture of h as a func- 
tion of horizontal position x and time t. We do this first, and later extend our discussion 
to determine shape and spreading kinetics in the (2 + l)-dimensional case. As illustrated 
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in Fi g. [H the z dependence of p z is determined by Frank elasticity ( DeGennes fe Prostl 



1993) with bo undary conditions p z = at z = an d p z = d x h <C 1 at the free surface 
z = /i, so that ( Sankararaman fc Ramaswamy||20"09l ) 



Pz = (z/h)d x h. (2.4) 

The z component of the equation of motion (|2.3[) reads ~o~od x p z = d z P where P is the 
pressure. Expressing p z in terms of h via (|2.4|) . using the Laplace pressure boundary con- 



dition with surface tension 7, and solving for P we find (jSankararaman fc Ramaswamv 
2009) 



P = P Q - jdlh + P a (z) (2.5) 

with P a {z) = — (ao/2)(z 2 — h 2 )(h" /h — h' 2 /h 2 ) where a prime (') denotes d x . In the 
following we ignore this active contribution to the pressure, which within the lubrication 
approximation gives only a negligible contribution to the equation of motion below. 

The equation of motion is obtained by integration of the x component of Eq. (12.3[) . 
using (|2.5j) . the active stress from (|2.1j) and including in the stress the contribution of 
nematic elasticity with the orientation profile (|2.4[) and a single Frank constant K. 

^Ul = y (idlh - dj-^ - ^d x h) (2.6) 

where uj_ is the inplane velocity averaged over the thickness of the film and F e i = 
(K/2) J dx(d x h) 2 jh is the nematic elastic free energy. Combining Eq 12.61 with the mass 
conservation Eq. 12.21 one obtains an independent equation for the height field h. 

d t h + — V • ( 7 /i 3 VV 2 /i - V/^V— ^ - a h 2 Vh) = (2.7) 
3/U oh 

where the gradient is taken here is the x direction. 

The effects of activity on spreading enter at the same order in gradients as those of 
gravity, but with a different dependence on height, essentially because of the existence 
of an intrinsic stress scale do- Indeed, the contribution of activity to (|2.6|) defines an 
active velocity v ac t — —Mhd x P actl with the hydrodynamic mobility Mh — h 2 /(3r]) and 
an active disjoining pressure 

Pact = 0-0 In h/ h Q (2.8) 

where ho is an arbitrary reference height. 

The shape of the static drop configuration that results at the end of spreading is 
particularly easy to understand in the 1 + 1-dimensional picture. In the static limit f|2 . 6|) 
integrated once with respect to x yields 

-ydlh - o-q In - ^ = const. = (2.9) 
ho oh 

where ho is the asymptotic maximum height for an infinite drop. The integration constant 
on the right hand side, i.e., the negative of the pressure jump at the top of the drop, 
has been taken to be zero, assuming a macroscopic flat drop. For a finite size drop, this 
pressure is obtained from the volume conservation constraint. Multiplying by d x h and 
integrating once again with respect to x yields, after some algebra, 

J(X + 7)(d^) 2 " Mhhi ^--h) = const. = Iq 2 (2.10) 

where 8 is the Young contact angle that the drop would achieve in the absence of activity. 
It can be seen as the effective contact angle made by the drop at the solid surface, 
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Figure 3. Height profile of a spread extensile drop, flat on top with a precursor region of 
quadratic variation at the limbs, obtained from a numerical solution of (|2 . 9|) for drop of finite 
volume fixed by a pressure of .0001 in units of ao- 



for heights approaching the value hx ~ K/j below which nematic elasticity becomes 
important. Above this height we can neglect K/h, and (|2.10|) takes the form of energy 
conservation for a particle of mass 7, position h(x) as a function of time x, moving in a 
potential 

V{h) = -a (h]n^- -h). (2.11) 
ho 

which vanishes at and has a maximum at ho- The steep growth of V(h) from h = 
implies a very slow progress of the effective particle, which means the height changes 
very slowly with x for h > hx, i.e., a flat free surface. The maximum height corresponds 
to the classical turning point for motion in the potential V(h) (|2 . 1 1 1) . d x h = at h = ho- 
From (|2.10l) this means 

7 6> 2 

ho = ~. (2.12) 

The drop is flat with a thickness ho over most of its area. At its edge, the thickness 
decays over a thickness t ~ ho/0 ~ j8/ao- The condition for the drop to be flat when 
it reaches equilibrium is that its final radius R be much larger than the size £ of the 
edge. For h — > the shape of the spreading drop is driven by the Frank elasticity of the 
polarization in (|2.10|) : 

(2-13) 

so that h ~ (jd 2 /AK)x 2 , growing quadratically from zero near the left edge which we 
have take n to be at 0, and similarly for the right edge. The resulting height profile is 
sketched in Fig. [3] 

The profile calculated and presented in (|2.11l) - (|2.13l) and Fig. [3] is consistent only if 
the height of the drop ho is larger than the thickness hx = if/7 below which nematic 
elasticity becomes relevant. This holds true if Q 2 > aoK/j 2 . If this condition is not 
satisfied the wetting properties of the drop are not dominated by the active stress but 
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Figure 4. Height profile of a spread contractile drop, higher and with a smaller base than in the 
extensile case (Fig. [3]), obtained from a numerical solution of (|2.9[) with the sign of the active 
stress reversed. 



by the nem atic elasticity. This r egime has been studied both experimentally and theo- 
retically bv lPoulard et al. (2006), and we do not consider it any further. As mentioned 
in the introduction, a thin film of active fluid becomes u nstable via a Freedericksz-like 
instability and flows spontaneously ( Voituriez et al. 2005) if its thickness is larger than a 
critical thickness h c ~ (K/o-q) 1 / 2 . A non- flowing steady state of a spread active drop can 
therefore exist only if its height ho < h c . To observe the effects we are discussing we must 
then have \ik < h c , i.e., 9 A < o-qK/^ 2 . Note that this condition is compatible with the 
condition that the wetting is dominated by the active stresses only if o-qK/^ 2 <C 1. We 
assume that this condition is satisfied in the following so that that the thin film formed 
by the active drop does not show any Freedericksz instability. 

For a contractile drop, ao < 0, the situation is less illuminating. For h < hx the height 
profile is as in the extensile case, viz., quadratic growth of the height, as the active stress 
do plays no role and the Frank energy dominates. Thereafter, we repeat the analysis of 
motion in the effective potential V(h) (|2.11j) . which now decreases with an infinite slope 
at h = 0, passes through a minimum at ho and increases thereafter. The height profile is 
thus higher and fatter than that of a sessile drop in the absence of activity. 

We now consider the spreading kinetics of a 1 + 1-dimensional active drop in its late 
stages where its thickness is smaller than the critical thickness h c for the appearance of 
spontaneous flow. The equation of motion for a spreading drop is given by Eq l2.7l where 
we neglect the nematic elasticity term that only plays a role at very small thicknesses. 
This spreading equation can be discussed at the scaling level. Except for the very edge 
of the drop, the surface tension term is small compared to the active stress term. This 
leads to dho/dt ~ —aoh-o/R 2 where ho is the maximum thickness at the center of the 
drop and R its radius. The two-dimensional volume of the drop Q ~ hoR is conserved. 
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We therefore obtain the spreading law 

i?M~(^) 1/4 (2.14) 

2.2. Three-dimensional treatment: topological defects 

We now consider the real three dimensional geometry and discuss the spreading kinetics 
of an active drop. The anchoring conditions that we impose can be satisfie d only if there 



exist defects in the nematic polarization field (jDeGennes fc Prostl 119931 ). We consider 
three types of defects: an aster-like cylindrical wedge disclination located at the center 
of the drop with polarization pointing radially from the defect line; a vortex-like wedge 
disclination located at the centre of the drop with polarization aligned in concentric 
circles in the azimuthal direction around the defect line; and two point-like half-defects 
located on the substrate on opposite sides of the drop as shown in Fig. [2l The first 
two topologies are simpler because the active drop remains axially symmetric but it is 
now well established that aster-like line defects are unstable and degenerate into two 
point-like half-defects. If the two defects remain on the vertical axis, the deformation of 
the polarization is virtually the same at that obtained for the original line defect. It is 
more likely that the two point-like defects migrate to opposite sides of the drop on the 
substrate which is our last topology. 

For the aster-like line defect (Fig. [5^), the polarization is mostly radial with a small 
component along the z direction. This is very similar to the 1 + 1 geometry studied 
above and if we ignore the core of the defect, the derivation of the equations for the 
local thickness of the drop exactly follows the same lines. The final shape of the drop 
is a circular pancake with a thickness given by equation 12.121 The equation governing 
the thickness during spreading is still Eq. 12.71 and a scaling analysis imposing a constant 
volume of the drop fi ~ hR 2 leads to the spreading law: 

««)~(^)" 8 (2.15, 

For a vortex- like line-defect (Fig. [2b), the polarization is always strictly azimuthal. In 
polar coordinates (r, 9), it is always along the 9 direction and the only non vanishing 
component of the active stress is a e a B = — ctq. The equation of motion in this case written 
in polar coordinates is 

i 2 

fiu ± = — ( 7 VV 2 H- r) = (2.16) 

3 r 

where f is the unit vector along r. The active stress appears here as similar to a centrifugal 
force with an unusual r dependence ~ 1/r instead of ~ r. 

The final stationary shape of the drop can be calculated explicitly by imposing a vanishing 
velocity on (|2.16l) : 

h = -^r 2 logr/R - ^(i? 2 - r 2 ) + A (i?2 _ ^ (2 . 17 ) 

where R is the radius of the drop. This profile is shown in Fig. \5\ The thickness has 
a minimum at the center of the drop where the polarization gradient is large. It has 
a maximum at an intermediate radius r = i?exp — (2j9/aoR). There does not seem to 
exist any simple scaling regime for the spreading kinetics of the drop in this case except 
at early times when the the active effects are negligible and the drop spreads as a simple 



A Drop of Active Matter 9 

h 



I \ 0.15 




/ \ 10 




/ o.oV 





-10 -5 5 10 



Figure 5. The final shape, from Eq. (|2.17(l . of an active drop in which the director field is in 
the vortex configuration, oriented strictly in the azimuthal direction, with a singularity on a 
vertical line segment at the centre of the drop. 



liquid drop R(t) ~ i 1 / 10 . The study of the spreading kinetics would require a numerical 
study of the spreading equations that we postpone for future work. 

The third topology with opposite defects on the two sides of the drop, as in Fig. Etc), is 
more complex since in this case the drop is anisotropic and spreads preferentially in the 
direction of the two defects and its radius in this direction R x is larger than its radius 
in the perpendicular direction along the planar substrate R y . As above, we suppose 
here that the polarization field follows adiabatically the deformation of the drop and 
reaches rapidly a local equilibrium where the molecular field conjugate to the polarization 
vanishes. When the anisotropy R x /R y is large, the polarization can be estimated from 
nematic elasticity as p x ~ 1, p y ~ —xy/R 2 and p z = {z/h)d x h over most of the volume 
of the drop. Inserting these results in the expression of the active stress one derives a 
dynamic equation equivalent to (|2.7p . Estimating the order of magnitude each term in 
this dynamic equation, we obtain two scaling relations (dR x /dt)/R x ~ (dR y /dt)/R y ~ 
aoh 2 /R 2 . Together with volume conservation hR 2 — constant, these equations do not 
suffice to obtain a growth law for the spreading drop as they do not give information on 
the anisotropy: this would require the knowledge of the prefactors of the scaling relations. 
Two possible limiting cases are (i) finite, bounded anisotropy, where the active drop is 
characterized by a single length-scale given by (|2.15|) and (ii) infinite anisotropy, where 
the drop can be considered as one dimensional with a spreading law given by (|2.14|) . 
Further progress requires a numerical solution, which we do not attempt here. 
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3. Summary and outlook 



In closing, we summarise our approach and results, and offer a perspective on future 
directions. From the point of view of fluid mechanics, the distinguishing feature of sus- 
pensions of self-driven particles is the permanent uniaxial stress with which each particle 
is endowed. Long-range orientational order organises the individual active stresses, trans- 
mitting their effect to macroscopic scales and leading to consequences for large-scale fluid 
flow. The spreading of a liquid drop on a solid surface, driven by gradients in the dis- 
joining pressure, is a fundamental problem in fluid mechanics. Active stresses provide 
a hitherto unexplored contribution to disjoining pressure. In an active suspension with 
nematic liquid-crystalline order, the deviatoric active stress is determined by the local 
ncmatic director which, in our treatment, is linked to the shape of the drop through the 
requirement that there is no normal component of the orientation field at any bounding 
surface with which the suspension is in contact. Thus at the base of the drop the director 
is purely horizontal, while variations in the free-surface height tilt the orientation. The 
resulting inhomogeneities in the active stress generate flows that drive the spreading pro- 
cess with novel growth laws and drop shapes. The closed geometry of a drop constrains 
the director field, inducing topological defects. The interplay of spreading and defects 
has not been studied even in conventional equilibrium nematics, as far as we know. We 
have studied 1 + 1 and 2 + 1 dimensional geometries. In a 1 + 1-dimensional geometry, 
corresponding to a distended drop, we are able to predict an anomalously flat final shape 
and a t 1 ^ 4 spreading law for extensile active stresses. For a full 2 + 1-dimensional case, 
scaling alone offers only partial information for the case where a pair of half-defects are 
situated at the ends of the drop. For the case of a single defect line at centre, we find a 
^l/e spreading law for a radial "aster" , and arrested spreading with a height profile with 
a non-monotone dependence for a vortex configuration, in which the director is strictly 
azimuthal. The complete 2 + 1-dimensional spreading law requires a numerical solution 
of the active thin-film equations. 



Before we end, let us note some of the simplifying assumptions we have made, some 
of which need to be relaxed in a more complete treatment. First is the approximation 
of a small contact angle. For our boundary condition, namely, filaments parallel to the 
bounding surface with which they are in contact, the effects of activity on spreading are 
strongest in the limit of a small contact angle. A contact angle close to 90° would mean 
that filaments lying along the free surface but close to the base would push vertically 
rather than laterally, suppressing the contribution of active stresses. Elastic anisotropy, 
which we have ignored, should le ad to spiral defect configurations which should rotate 
spontaneously ([Kruse et all 120041 ) . If we take into account the polar nature of the order 
parameter, even a pure vortex should rotate, and a macroscopically polar active drop 
should translate. We have not examined so far the stability of active drops with respect to 
azimuthally varying perturbations. Finally, in vivo realizations of active drops in general 
have complicating features that will need to be suppressed experimentally or, eventually, 
included in a comprehensive theor etical analysis. For ex ample, cell division will modify 
spreading in bacterial suspensions (jBees et aZ.ll2000l . 120021 ). From a long-term perspective, 
our work is of relevance to an understanding of the s preading of a single cell. However, 
in that setting, treadmilling ()Verkhovskv et al. Il999l) . that is, the polymerization and 
depolymerization of actin filaments, will play an important role. Notwithstanding these 
limitations, we look forward to experimental tests of our prediction s, on microorganism 
suspensions, extracts, or artificial self-pr opelled systems (see, e.g. 
Howse et al\ 120071 : iBartolo k Laugall2010l ). 



Paxton et all [2004; 
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